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Abstract 

We test the concepts of renormalized charge and potential satura- 
tion, introduced within the framework of highly asymmetric Coulomb 
mixtures, on exactly solvable Coulomb models. The object of study is 
the average electrostatic potential induced by a unique "guest" charge 
immersed in a classical electrolyte, the whole system being in ther- 
mal equilibrium at some inverse temperature (3. The guest charge is 
considered to be either an infinite hard wall carrying a uniform sur- 
face charge or a charged colloidal particle. The systems are treated 
as two-dimensional; the electrolyte is modelled by a symmetric two- 
component plasma of point-like ite charges with logarithmic Coulomb 
interactions. Two cases are solved exactly: the Debye-Hiickel limit 
/?e^ — > and the Thirring free-fermion point j3e^ = 2. The results at 
the free-fermion point can be summarized as follows: (i) The induced 
electrostatic potential exhibits the asymptotic behavior, at large dis- 
tances from the guest charge, whose form is different from that ob- 
tained in the Debye-Hiickel (linear Poisson-Boltzmann) theory. This 
means that the concept of renormalized charge, developed within the 
nonlinear Poisson-Boltzmann theory to describe the screening effect of 
the electrolyte cloud, fails at the free-fermion point, (ii) In the limit of 
an infinite bare charge, the induced electrostatic potential saturates 
at a finite value in every point of the electrolyte region. This fact 
confirms the previously proposed hypothesis of potential saturation. 
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1 Introduction 



Asymmetric classical Coulomb mixtures, such as highly charged colloidal 
or polyelectrolyte suspensions, in the strong coupling regime, have attracted 
much attention in the last years (for a review of phenomenological approaches 
built on the base of mean- field theories, see ref. PJ). The concept of renor- 
malized charge has been introduced within the Wigner-Seitz cell models 
to describe an effective interaction between highly- charged "macro-ions" as 
a result of their strong positional correlations with the oppositely charged 
"micro-ions" 013 1313 El- 

The concept of renormalized charge can be documented in the infinite 
dilution limit of colloids [3 |Hl 13 113 HI] ■ The simplified model consists of 
a unique colloidal particle idealized as a hard sphere of radius a carrying 
charge Ze, Z is the valence and e the elementary charge, at its centre (any 
spherically symmetric charge distribution inside the colloid can be repre- 
sented in this way). The colloid is immersed in an electrolyte modelled by 
a symmetric two-component plasma (TCP) of elementary ±e (if possible, 
point-like) charges. The system is defined in an infinite i/-dimensional space 
of points Y E W , having for simplicity vacuum dielectric constant e = 1. 
The interaction energy of two particles of charges q and q' at the respective 
spatial positions r and r' is given by gg'0(|r — r'|), where 0, the Coulomb 
potential induced by a unit charge, is the solution of the Poisson equation 

A0(r) = -s,5{v) (1.1) 

Sy is the surface area of the i/-dimensional unit sphere. This definition of the 
i/-dimensional Coulomb potential maintains many generic properties (e.g., 
screening sum rules ^2|) of "real" three-dimensional (3D) Coulomb systems 
with the interaction potential 0(r) = 1/r, r G R^. In particular, in 2D, 

0(r) = -ln(|r|/ro), r G (1.2) 

where tq is a free length scale. Thermal equilibrium is treated in the grand 
canonical ensemble characterized by the inverse temperature (3 = 1/ [kT) and 
by the couple of bulk particle fugacities for electrolyte ±e-charged particles, 

= = z. The corresponding average bulk densities will be denoted 
by 72+ = n_ = 72/2, n is the total number density. For the classical case 
of point-like particles, the singularity of the Coulomb potential 0(r) at the 
origin r = often prevents the thermodynamic stability against the collapse 
of positive-negative pairs of charges: in 2D, for small enough temperatures; 
in 3D, for any finite temperature. 

We shall explain the concept of renormalized charge on the 3D mean- 
field theories, valid in the high-temperature limit and free from the collapse 
problem of point-like charges. Let us fix the colloidal particle at the origin 



2 



and denote by ipi^) the induced average electrostatic potential at point r. 
Inside the colloidal hard-core region < r < a, iplr) satisfies the Poisson 
equation 

Aip^ir) = -A7rZe5{r) (1.3) 
Inside the electrolyte region r > a, the Poisson equation takes the form 

A^>(r) = -47rp(r) (1.4) 

where 

p(r) = e[n+(r) -n_(r)] (1.5) 

is the charge density of electrolyte particles. Within the ordinary mean-field 
approach, the average particle densities at a given point are approximated by 
replacing the potential of mean force by the average electrostatic potential 
at that point, n±(r) = n±exp[^Peip^{r)]. Eq. p.4|l then reduces to the 
non-linear Poisson-Boltzmann (PB) equation 

A'?/^>(r) = 47rensinh [f3eip-^{r)] (1-6) 

Debye and Hiickel proposed a linearization of this equation by consider- 
ing the small-argument expansion of smh.{f3eipy) ~ f3eip^ valid in the high- 
temperature limit. The linear PB equation then reads 

(A-«:2)^>(r) = (1.7) 

where k = ^J4:^T(3e^n is the 3D inverse Debye length of the electrolyte par- 
ticles. Due to the spherical symmetry of the problem, A = 9^ + {2/r)dr- 
Eqs. (jl.3j) and (jl.7|) . subject to the condition of regularity at r — > cxd and 
the boundary conditions of continuity of the electrostatic potential and the 
electric field across the colloid surface r = a, imply that, in the linear Debye- 
Hiickel theory, 

Mr) = Ze(---^) (1.8) 
Vr 1 + K.aJ 

^>(r) = exp[-.(r-a)] 

The charge density of electrolyte particles is given by 

p(r) = —Pe^nipyir), r > a (1-10) 

Strictly speaking, the linearization of Eq. ()1.6p is only valid if l3etp{r) << 1, 
and this is indeed true for asymptotically large distances r where the screened 
potential vanishes. As a consequence, at large r, the solution of the non- 
linear PB equation (jl.6j) must also take the Yukawa form 

, , , exp \—K,(r — a)] , ^ 

V^>(r)~A— ^-^ ^ ^, r^oo (1.11) 
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with A being an r-independent constant. The non-hnearity of the problem is 
reflected only through an effective boundary condition at the colloid surface 
determining the constant A. Comparing with ()1.9|) one sees that A is 

related to the renormalized colloid valence, Zj-en, as follows 



The renormalized charge reflects the screening effect of the electrolyte cloud, 
and can further be used to establish an effective interaction for a system 
of guest charges immersed in the electrolyte. Since the exact solution of 
the non-linear PB equation for the sphere geometry is not available, 
the constant A can be determined only approximately, i.e., in the limiting 
case Ka » 1, by matching with the exact solution of the non- linear PB 
for the charged-plane geometry E]. An important feature is that, as 
expected from the Manning condensation theory [2^ , the renormalized valence 
saturates at some finite value Z^^* when the colloidal bare valence Z goes to 
infinity. 

More refined approaches, which go beyond the mean-field approximation 
and incorporate electrostatic correlations among the electrolyte particles, 
were developed in refs. ^151 114j. Monte-Carlo simulations ^H] indicate the 
existence of a maximum in the plot of the renormalized charge versus the 
bare colloidal charge. 

As was correctly mentioned by Tellez and Trizac ^^1; the definition of 
renormalized charge requires that the average electrostatic potential behaves 
far from the colloid as it would within the linearized Debye-Hiickel theory, up 
to the constant prefactor. This is not at all ensured for a finite temperature. 
To avoid this artificial limitation in the saturation problem, one considers 
the possibility of a more general phenomenon of potential saturation: the 
question is whether, in the limit of the infinite bare colloidal charge Ze ^ oo, 
the induced average electrostatic potential tp^^^{r) is finite inside the whole 
electrolyte region r > a ^. We would like to emphasize that the potential 
saturation, if it exists, is the pure non-linearity effect: there is no potential 
saturation within the linear Debye-Hiickel theory, see Eq. (II. 9|) . 

It is well-known that the linearized Debye-Hiickel theory correctly de- 
scribes the small-coupling (high-temperature) limit /3e^ — > in the sense 
that the basic screening properties of the charged system in the conducting 
regime are preserved. In the present case one can readily verify by using 
Eqs. (jl.9|) and (|1.1(J|) that the screening cloud of the electrolyte particles 
compensates exactly the bare charge of the "guest" colloid: 



based on simple electrostatics, the potential is infinite at the colloid surface r — a 



A 



1 + Ka 
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This sum rule no longer holds within the non-linear PB theory. This serious 
deficiency of the non-linear PB theory is "tolerated" in various phenomeno- 
logical approaches because of the predicted charge renormalization which 
is in a relatively good agreement with available experimental data. 

The present work aims to put the concept of charge renormalization and 
the hypothesis of the electric potential saturation on a rigorous basis. As 
a test model for the electrolyte, we use the 2D symmetric TCP of point- 
like ±e charges with logarithmic pairwise interactions ()1.2|) . The 2D plasma 
of point-like charges is stable against the collapse of positive-negative pairs 
of charges provided that the corresponding Boltzmann factor r~^^ is inte- 
grable at short distances in 2D, i.e. for the (dimensionless) coupling constant 
/3e^ < 2. In this stability range of couplings, the equilibrium statistical me- 
chanics of the plasma (the bulk thermodynamics, special cases of the surface 
thermodynamics and the large- distance behavior of the two-body correla- 
tion function) is exactly solvable via an equivalence with the integrable 2D 
Euclidean sine-Gordon field theory (for a short review, see ref. ^|). The 
complete exact information about correlation functions is available in two 
special cases: in the high-temperature Debye-Hiickel limit /5e^ — > 0, and just 
at the collapse point /?e^ = 2 |T71 1^ which corresponds to the free-fermion 
point of an equivalent 2D Thirring model (although the free energy and the 
particle density diverge at the collapse point, the truncated Ursell correlation 
functions are finite). We examine the above-outlined problems in these two 
exactly solvable cases, for two particular geometries: the charged line and 
the charged circular colloid. Based on the exact results we show that the 
concept of renormalized charge does not apply to the studied 2D microscopic 
Coulomb system. On the other hand, the anticipated phenomenon of electric 
potential saturation is confirmed. 

The paper is organized as follows. Section 2 deals with the charged-line 
geometry. Section 2.1 is devoted to a short recapitulation of the Debye- 
Hiickel limit. In Section 2.2, the known exact results at the free-fermion 
point jlHI are analyzed from the point of view of the studied subjects. Section 
3 deals with the colloidal- charge geometry. As before. Section 3.1 concerns 
the Debye-Hiickel limit, while Section 3.2 is devoted to an original derivation 
of the exact solution at the free-fermion point. A recapitulation and some 
concluding remarks are given in Section 4. 

2 Charged line 

We consider an infinite 2D space of points r E B? defined by Cartesian co- 
ordinates (a;, y). The half-space a; < 0, impenetrable to particles, is assumed 
to be a vacuum hard wall. The electrolyte of ±e point-like charges is con- 
fined to the complementary half-space a; > 0. The model interface is the 
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line localized at x = 0, along the y-axis. The line, which carries a uniform 
charge ere per unit length, models an electrode. There is another electrode 
of opposite charge density localized at a; = +00. The electrostatic potential 
induced by the two electrodes is for a; < and —2'7iaex for x > 0. The 
boundary condition for the electric field reads 



dx 



(2.1) 



2.1 Debye-Hiickel limit 

The average electrostatic potential at distance x form the interface satisfies 
the 2D Poisson equation 

-fV^ = -27rp(x), X > (2.2) 
dx^ 

In the spirit of the Debye-Hiickel theory valid in the limit /3e^ 0, the 
charge density of the electrolyte particles is approximated, in close analogy 
with Eq. (jLlUp . as p{x) ~ —Pe^nip{x). The linear PB equation thus reads 

^{x) = 0, X > (2.3) 



dx^ 



where k = y/2TT(3e^n is the 2D inverse Debye length. The solution of ()2.3|) . 
subject to the requirement of regularity at x — > 00 and the boundary condi- 
tion (j2.ip . takes the form 

27rcre 

ipix) = exp(— Kx), X > (2.4) 

K 

The consequent charge density 

p(x) = — acK exp(— Kx), x > (2.5) 
fulfills the following analogue of the screening sum rule ()1.13j) : 

PCX) 

ae+ / dxp(x) = (2.6) 



2.2 The free-fermion point 

Since the particle density diverges at the collapse point /5e^ = 2, the avail- 
able thermodynamic parameter is the fugacity z. It will be considered in a 
rescaled form, m = 2nroz [tq is the length scale considered in p.2j) ]. having 
dimension of an inverse length. The density profiles of electrolyte particles 
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near the charged plain hard wall were obtained in ref. ^Hl by solving the 
Green-function problem of the corresponding boundary Thirring model at 
its free-fermion point. The result is 

2 

771 

n±{x) = 7l± + - — [K2{277LX) — Ko{277lx)] 



2 

m 



27r 

2 

771 
2^ 



1 

+ 



2777X ' {277^xy\ ^-P(-2-^) (2.7) 

/■'^^'"^ dt / ^ N 

/ , = exp — 2v^ +t a; 

Jo + t^-t ^ ^ 

where Ki are the modified Bessel functions of order /. The divergent bulk 
particle densities n+ = n_ = n/2 can be regularized, e.g., by considering 
a small hard core around each particle ^Hj; since we are interested in the 
charge density p{x) defined by the difference e[n+(x) —7i_{x)], we avoid this 
regularization procedure. After some simple algebra, one gets 

p{x) = / dt Vm2 + ^2 exp (-2^771^ + f^x) (2.8) 

TT JO ^ ^ 

It is easy to check that the screening sum rule ()2.6|) is fulfilled by this charge 
density. 

The corresponding electrostatic potential, determined by the Poisson 
equation ()2.2|1 and the requirement of regularity at x — ^ oo, reads 

e /■27r(T 



tPix) = - / ^== exp (-2Vm^+¥x) (2.9) 

The boundary condition ()2.1|) is evidently satisfied for this potential. In order 
to obtain the large-x expansion of ipix), we first make in ()2.9|) a change of 
the integration variable t into u = x[^l + (t/m)2 — 1], and then expand the 
integrated function in powers of 1/x. In the leading order, 

■^(x) ~ - (77-7 ) exp(— 2ma;), x —>■ (yo (2.10) 



TTIX 



for all 0" 7^ [cr = implies the trivial result tplx) = 0]. The independence 
of the leading asymptotic term ()2.10|1 on the (nonzero) charge density a is 
a special feature of the present geometry. Comparing ()2.1()j) to the Debye- 
Hiickel result ()2.4j) characterized by the pure exponential decay in x, and 
identifying the respective inverse lengths 2m and k, one sees that the large-x 
behaviors differ one from the other. The idea of renormalized charge density 
thus fails. On the other side, considering the limit of the dimensionless ratio 
cr/m ^ 00 in fl2.9j) . one has explicitly 

^Ij^^\x) = |ifo(2mx) (2.11) 

It follows from the basic properties of Kq that < ip^^^lx) < 00 for all x > 0, 
in full agreement with the saturation hypothesis. 
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3 Charged circular colloid 



As in the 3D case discussed in the Introduction, we fix at the origin one 
colloid of charge Ze and disk hard core with radius a. There is an infinite 
2D TCP of ±e point-like charges in the complementary outer space. The 
analogue of the boundary condition ()2.1|1 for the electric field now reads 



dr 



Ze 



(3.1) 



3.1 Debye-Hiickel limit 

Inside the colloidal hard-core region < r < a, the electrostatic potential 
satisfies the 2D Poisson equation 

A^/'<(r) = -27rZe(5(r) (3.2) 

Inside the electrolyte region r > a, the consideration of the linear mean-field 
prescription p(r) ~ — /3e^n^/'(r) in the Poisson equation implies 

(A-«:2)^>(r) = (3.3) 

Due to the circular symmetry of the problem, A = 5^ + {l/r)dr- Eqs. ()3.2p 
and ()3.3|) . subject to the requirement of regularity at r — > oo and the usual 
boundary conditions of continuity across the colloid boundary r = a, imply 



KaKAKa) 



(3.4) 



^>{r) = . M^r) (3.5) 

KaKi[Ka) 

The boundary condition ()3.H) is trivially satisfied. Note that, after defining 
the surface charge density ae = Ze/{2Tra) and going to the limits a,r ^ oo 
with a fixed difference r — a = x > 0, ()3.5|) reduces to the straight-line result 
()2.4|1 as it should be. At large r, using the asymptotic formula for Kq fi9\ . 
the average electrostatic potential ()3.5j) behaves like 

yjy{r) — TT-, — r ( 7^ — ) exp(— /tr), r — > oo (3.6) 



KaKi{Ka) \2Kr 
The electrolyte charge density, given by 

P^'^^ = ~n JF~( — V^o(«'^), r>a (3.7) 

fulfills the screening sum rule 

roo 

Ze+ dr27rrp(r) = (3.8) 

Ja 
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3.2 The free-fermion point 

According to the general formalism established in ref. fTli I18j . in order to 
obtain density profiles of electrolyte ±e particles at coupling /3e^ = 2, one 
has to solve the Green function problem of a 2 x 2 matrix G(r, r'). Its matrix 
elements Ggg'(r, r') (g, g' = ± denote the charge sign) are determined by a 
system of four coupled partial differential equations (PDE), written in a 2 x 2 
matrix notation as follows 



1 + _ 1 - 



a^dx + (j'^dy + m+(r) h m_(r 



2 ' ' 2 



G(r,r') = I5(r-r') (3.9) 



Here, 1 and a* (i = 1, 2, 3) denote the 2x2 unit and Pauli matrices, respec- 
tively, and 

mg(r) = m(r) exp [— 2gt'(r)] , g = ± (3.10) 

is the position-dependent (rescaled) fugacity for some external electric po- 
tential v{r) (in units of e); a nonelectric potential which acts in the same 
way on both kinds of particles, like an impenetrable hard wall or core, is 
described by the region- dependence of m(r). Four Eqs. (j3.9|) split into two 
independent sets of equations, the one for the pair (G++, and the other 

for the pair {G , G'+_). We shall present a detailed derivation of the results 

for the pair (6*++, G^+), given by 

m+(r)G++(r,r') + (9,-i9,)G_+(r,r') = 5(r - r') (3.11) 
(9,. + ia,)G'++(r,r')+m_(r)G'_+(r,r') = (3.12) 

Based on a similar derivation procedure, we shall only present the final results 

for the pair (G ,G'_|__). As concerns the boundary conditions, since Eq. 

(|3.9|) is a first-order system, all matrix elements Ggg'(r, r') must be continuous 
when crossing a boundary between two different regions. The requirement 
of regularity is obvious. The one-particle densities are given by 

ng{r) = mg{r) lim G,,(r, r'), q = ± (3.13) 



In the bulk regime with m(r) = m for all points r E R , one has 

TTl 

GUr,r') = —Ko{m\r-r'\), q = ± (3.14) 
Zn 

so the one-particle densities diverge logarithmically as r' — > r in ()3.13j) . This 
divergence can be suppressed by introducing a short-distance cut-off R, 



n+ = hm — Kn{mR) ~ — 



' mi?) 



(3.15) 



where G is Euler's constant. 
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In the case of the studied model, the colloidal particle at the origin induces 
the electrostatic potential (in units of e) v{r) = — Zln(r/ro). The rescaled 
fugacity m(r) = inside the hard-core region < r < a and m(r) = m in 
the electrolyte region r > a. Thus, 

{0 for < r < a 

(3.16) 
m(r/ro)^^^ for r > a 

For our purpose it will be sufficient to consider the source point r' in Eqs. 
()3.11|) and ()3.12p to be localized only in the electrolyte region, so, without 
writing it explicitly, in what follows we shall assume that r' > a. As concerns 
the point r, let us first assume its localization in the colloidal hard-core 
region, i.e. r < a. Taking mj-{r) = in Eqs. (|3.11|) and (|3.12j) . one gets 

id, + idy)G++ir,r') = 0, r<a (3.17) 
(9,-i9,)G'_+(r,r') = 0, r<a (3.18) 

This means that, as functions of r, 6*++ depends only on z = rexp{i(p) and 
G_4. depends only on z = r exp{—i(f). The general solutions of Eqs. ()3.17p 
and ()3.18|) . regular at r = 0, can be therefore written as the expansions of 
the forms 

m °° 

G'++(r,r') = — ^/i(mr',(y9')(mr)'exp(i/(^), r<a (3.19) 
2^ 1=0 

m °° 

G'_+(r,r') = — V/ii(^^'.V^')(^0'exp(-i/(^), r<a (3.20) 

where the functions {fi,hi} are determined by the boundary conditions at 
r = a. When the point r is localized in the electrolyte region, i.e. r > a, 

taking m±(r) = m(r/ro)^^^, we first express by using Eq. ()3.12j) G ^ in 

terms of 

G_+(r, r') = -- (-Y^ (9, + idy) G++(r, r'), r>a (3.21) 

m xtq/ 

The consequent substitution of into p.llj) implies the only PDE de- 
termining After lengthy but simple algebra, in terms of the auxiliary 
two-point function 

^?++(r,r') = -- (-)^G++(r,r') (-] (3.22) 

this PDE is obtained in the form 

(-m2-//)^++(r,r') =5(r-r'), r>a (3.23) 
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where 

H = -A^- — {xdy - yd,) + — (3.24) 

It is clear that (7++(r, r') = (r|(— — H)^^\r') is nothing but the Green- 
function two-point matrix element associated with the one-particle quantum 
Hamiltonian H and the spectral parameter — m^. In polar coordinates (r, ip), 
the Hamiltonian ()3.24j) reads 

^ = -^-1^-1^-^^ + ^ (3.25) 

According to elementary quantum mechanics, the periodicity requirement 
under the shift ip ^ ip + 27t implies that the eigenf unctions of H have a 
trivial dependence on the angle ip: oc exp{il{p) where / = 0, ±1, ... is the 
"magnetic" quantum number. It follows from Eq. ()3.25j) that the radial part 
of the eigenfunction with a given / is then determined by 

1 d (l + Z)^ , , 

This is the radial Hamiltonian of a free quantum particle in 2D, where the 
presence of the colloidal charge Ze manifests itself as the shift of the quantum 
number / by the integer Z. The standard Green-function technique implies 
an explicit form of (?++. Using then the relation ()3.22|) . G'_|__|_(r, r') with 
r, r' > a is found to be 



G 



l=—oo 



X [Ii+z{mr<)Ki+z{mr>) + CiKi+z{'mr)Ki+z{mr')] (3.27) 

Here, /; and Ki are the modified Bessel functions, and r< (r>) is the smaller 
(the larger) of r and r' . is generated from 6*++ via Eq. ()3.2H) . For the 
special case a < r < r', it takes the form 

G_,(r,rO = E^exp[i(/+l)^-iV] 

X [-Ii+z+i{mr) + ciKi+z+i{mr)]Ki+z{mr') (3.28) 

The unknown constants {q} are determined by the requirements that G++ 
and G_+ be continuous at r = a. With regard to Eqs. ()3.19j) and ()3.2()|1 . 
one gets 

-Ii+zima) / Ki+zima) I < 
ci = { (3.29) 

Ii+z+i{'ma)/Ki+z+i{'ma) l>0 
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Applying the "summation theorem" for the modified Bessel functions [TU] 

oo 

Ko{m\r - r'\) = ^ exjp[il{ip - ip')]Ii{mr^)Ki{mr>) (3.30) 



(=— oo 

we conclude that, when both points r and r' are in the electrolyte region 
(r, r' > a), 



+ Eexp[i/(¥. - ¥.')]^±^i^T^^^+z(mr)ir,+z(mr') (3.31) 
-Eexp[i/(v.-^')]|^4^^/+z(mr)ir,+^(mr')| 

The matrix element G (r, r') can be deduced in the same way. It is 

expressible by formula 1)3.311) under the substitution Z —Z, which corre- 
sponds to the sign reversal of the colloidal charge. 

The densities of electrolyte particles at r > a follow from the relation 

^ Ii+i{ma) ^ Ii{ma) ^ 

n-i-{r) = n+ -\ > ; -\Kj{mr)\ > — ; A/ mr 

27r (ma) ^ 27r Ki {may 

(3.32) 

Using the Wronskian relation [T^ 

Ii{x)Ki+i{x) + Ii+i{x)Ki{x) = - (3.33) 

X 

together with the symmetry properties // = and Ki = K_i for integer /, 
the charge density in the electrolyte region is found in the form 

P{r) = E ^ f \r^f V ([Ki-i{mr)]^ + [Ki{mr)]^) , r>a 

2vra jz{ Ki_i[ma)Ki[ma) ^ ' 

(3.34) 

It stands to reason that p(r) = in the colloidal region r < a due to the 
hard-core repulsion. With the aid of the integral formula 

Arr\Ki{mr)f = y [Ki_,{ma)Ki+,ima) - [Ki{ma)]^) (3.35) 

and the recursion relation 

x[Ki^^ix) - Ki+,{x)] = -2lKi{x) (3.36) 
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taken a.t I = Z, the charge density (j3.34p can be shown to fulfill the screening 
sum rule (j3.8|) . 

The average electrostatic potential is given by the couple of Poisson equa- 
tions 

A^<(r) = -27rZe5{r), r<a (3.37) 
A^>(r) = -27rp(r), r>a (3.38) 

The circularly symmetric solution of Eq. (j3.37j) reads 



4'<:{r) = Ze 



—In ( - ) + const 



(3.39) 



The circularly symmetric potential of Eq. ()3.38|) fulfills the differential equa- 
tion 



d 
dr 



d^>(r) 
r- 



dr 



a Ki_i{ma)Ki{ma) 

(3.40) 

The first integration of Eq. ()3.4()|1 with respect to r can be performed easily. 
Using the regularity condition lim^^oo ^ dr^/'>(r) = 0, the integration formula 
of type ()3.35|) and finally the recursion relation ()3.36p . one arrives at 

dV>>(r) _ Cy. Ki^i{mr)Ki{mr) 
dr ^i^i Ki^i{ma)Ki{ma) 

Note the obvious fulfillment of the boundary condition ()3.1|1 . The subsequent 
integration of Eq. ()3.4H1 with respect to r is a bit more complicated problem. 
The regularity condition lim^^oo V'>(^) = has to be combined with the 
integral formula [derivable by using the relation ()3.52p ] 

i-i 



m dr' Ki^i{mr')Ki{mr') = {-l)^+^-Y,{-iy^i,j[Kj{mr)f, l>l 

Jr 2 .•_n 

(3.42) 



IX j is the Neumann factor: yUo = 1 and /ij = 2 for j > 1, to get 



= -^^Y.(-iyNfArna)\K,{mr)f (3.43) 
^ i=o 

where 

^ 1 1 
fAma) = V (-1)' — — r (3.44) 
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With regard to the Wronskian relation (|3.33|) . fj can be simphfied as follows 

fj{ma) -- 



z 

E(-i 

i=j+i 



Ii{ma) ^ Ii_i{ma) 



Ki{ma) Ki^iyma) 
z Izifna) 



Kz{ma) 

Thus, in the electrolyte region r > a, 



^_^ylj(jnal 
Kj {ma) 



(3.45) 



z-i 



V. ^ Ij{ma) 



[Kj{mr)f 



(3.46) 



Before analyzing the result ()3.4(i|l . let us recall some basic properties 
of the modified Bessel functions Ii{x) and Ki{x) {I = 0,1,...) when the 
argument x belongs to the interval < a; < oo. Ii{x) and Ki{x) satisfy the 
same differential equation. 



d'f ^ 1 df 



- 1 



dx^ X dx \ x^ ^ 
but exhibit different asymptotic behaviors: 



/ = 0, f = Ii{x) oiKiix) 



(3.47) 



Iiix) 



V2 



K,(x) 



TCX 



TT \l/2 



for X — s> oo 



and 



Iiix) 



1 fx 



l\ \2 



Ki{x) 



(/-I)! fx 



for X ^ 



(3.48) 



(3.49) 



except for the special / = case of -^'0(2^) ~ — ln(x/2) — C. Both Ii{x) and 
Ki{x) are positive for x > 0. In particular: Ii{x) starts from at x = 
[except for the special case of /o(0) = 1] and monotonously increases to 
infinity at x — 00; Ki{x) is infinite at x = and monotonously decreases to 
zero at X ^ 00. 

For finite Z, at large r, the average electrostatic potential ()3.46p behaves 

like 



Ij(ma) Izima) 
f^o^'Kj{ma) ^ Kz{ma) 



exp(— 2mr), r — >■ 00 

mr 



(3.50) 
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This asymptotic behavior differs from the large-distance prediction (|3.6|) of 
the Debye-Hiickel theory by the factor r~^/^, which is in contradiction with 
the concept of renormahzed charge. 

In order to resolve the saturation problem in the limit Z ^ oo, we first 
use the integral formula ()3.42j) to rewrite the last term on the rhs of Eq. 
flTinil as follows 

+e^^^^^m r dr' Kz-i{mr')Kz{mr') (3.51) 
Kz{ma) Jr 

The positivity of the modified Bessel functions then ensures that, for every 
Z, ip>{r) > in the whole region r > a. In order to establish an upper 
bound for ^/'>(r), we first use the relation jTHI 

Ki_,{x) + K,+^{x) = -2^^ (3.52) 
to establish the equahty 

POO poo 

m I dr' Kz-i{mr')Kz{Tnr) = [Kzimr)]^ — m / dr Kz{mr)Kz-\-i{'mr) 

Jr J r 

(3.53) 

The consideration of this equality in Eq. ()3.51|) leads to 

< ? E ^7^[^.(^0]^ r > a (3.54) 

In the considered region r > a, the inequality Kj{mr)/Kj{ma) < 1 implies 
that 

e ^ 

ip>{r) < - E Ij{fT^o)Kj{mr), r > a (3.55) 
2 j=-z 



In the limit Z ^ oo, the application of the summation theorem ()3.30|) with 
if = (f' finally gives 

0<^"'"\r) <^Ko{m{r-a)), r>a (3.56) 

The existence of the lower and upper bounds for tp^'"^^ in the electrolyte region 
confirms once again the validity of the potential saturation hypothesis ^3]. 



15 



4 Conclusion 



We have studied the average electrostatic potential induced by a unique 
"guest" charge immersed in an infinite electrolyte, the electrolyte being 
modelled by the classical two-component plasma of elementary ±e point- 
like charges. The primary motivation came from the predictions of the two 
basic 3D mean-field theories described in the Introduction: the Debye-Hiickel 
theory based on the linear PB equation and the non-linear PB theory. The 
important point is that both mean-field theories predict the same type be- 
havior of the induced electrostatic potential at asymptotically large distances 
from the guest charge, only the constant prefactors are different. Within the 
non-linear PB theory, this fact permits one to introduce the renormalized 
guest charge which involves the non-linear screening effect of the electrolyte 
cloud, and can further be used to establish an effective interaction for a 
system of guest charges immersed in the electrolyte. When the bare guest 
charge goes to infinity, the renormalized charge saturates at a finite value. 

In order to go beyond the mean-field methods, we have have tested the 
concept of charge renormalization on the 2D Coulomb models. These system 
have advantage of being completely solvable not only in the Debye-Hiickel 
high-temperature limit /3e^ — * 0, but also at a finite temperature, namely the 
Thirring free-fermion point jSe"^ = 2. Although just at this point the collapse 
of positive-negative pairs of point-like charges emerges, the charge-density 
profile in the electrolyte region (determining the average electrostatic po- 
tential through the Poisson equation) is a well-defined finite quantity which 
satisfies the electroneutrality sum rule. We have considered two geometries 
of the guest charge: the infinite hard wall carrying the uniform surface charge 
(Section 2) and the charged colloidal particle with a hard core (Section 3). 
For both geometries, the results at the free-fermion point can be summarized 
as follows. The asymptotic large-distance behavior of the induced electro- 
static potential differs from that predicted by the linear Debye-Hiickel theory, 
so the concept of renormalized charge does not apply. On the other hand, 
when the bare guest charge goes to infinity, the induced potential saturates at 
a finite value in each point of the electrolyte region. In the case of the infinite 
charged wall, the saturation potential was found explicitly, see Eq. ()2.11|) . In 
the more complicated case of the charged colloidal particle, lower and upper 
bounds for the saturated potential were derived, see formula ()3.56p . These 
results confirm that the potential saturation hypothesis jTH] is indeed true. 

It is an open question whether the failure of the concept of renormalized 
charge is restricted to the free-fermion point, or it takes place also in some 
subinterval of the stability region < /3e^ < 2. This question might be 
answered by exploring a form-factor approach to the integrable 2D sine- 
Gordon model, the field-theory equivalent of the 2D TCP. Without going into 
detail, the form-factor method serves as a tool for a systematic generation of 
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the large- distance asymptotic expansion for two-point correlation functions. 
Its recent developments seem to be applicable to the present problems. We 
shall proceed along this line. 
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